Cu NMR and NQR evidence for an unusual spin dynamics in PrCu2 below 100 K 
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We report the results of a "^^'^^Cu NMR/NQR study probing the intermetallic compound PrCu2. 
The previously claimed onset of magnetic order at 65 K, indicated in a /iSR study, is not confirmed. 
Based on our data we discuss different possible reasons for this apparent discrepancy, including 
a non negligible influence of the implanted muons on their environment. Competing dipolar and 
quadrupolar interactions lead to unusual features of the magnetic-ion/conduction-electron system, 
different from those of common intermetallics exhibiting structural or magnetic instabilities. 

PACS numbers: 71.70.Ej, 75. 20. En, 76.60.-k, 76.60.Es 



I. INTRODUCTION 

The Van Vleck paramagnet PrCu2 exhibits peculiar 
electronic and magnetic properties, whose complete un- 
derstanding is still to be achieved. The 4/ electron 
ground-state multiplet of the Pr'^+ ions is split into 
nine singlets by the crystal electric field. Nevertheless, 
an induced Jahn- Teller (JT) transition is observed at 
TjT = 7.6 K,— 1^ where the crystal structure changes from 
orthorhombic to slightly monoclinic^i^i^ An incommen- 
surate antiferromagnetic (AF) order below T/v = 54 mK 
among the Pr nuclear magnetic moments was claimed 
on the basis of magnetic susceptibility, specific heat, and 
neutron scattering measurements The respective mo- 
ments align along the a direction of the orthorhombic 
unit cell with a propagation vector ~q = (0.24,0,0.68) 
in the reciprocal lattice. Magnetization and suscepti- 
bility measurements indicate a large polarization of the 
4/ Pr magnetic-moments in the presence of an external 
field.— >^ At low temperatures the easy and hard magneti- 
zation axes are along the a and c direction, respectively. 
Quite surprisingly the two magnetization axes are ex- 
changed for fields exceeding 10 T, leading to a metam- 
agnetic transition.— The 4/ electron Pr-quadrupole mo- 
ments play an important role in the physics of PrCu2. 
In particular the JT transition is ascribed to a phonon- 
mediated interaction between the quadrupole moments,— 
whereas the metamagnetic transition has been ascribed 
to a change in the orientation of the spatial charge 
distribution!^ 

Based on the results of a muon spin resonance (/^SR) 
study, AF ordering of the locaHzed Pr 4/ electron mo- 
ments below 65 K has recently been suggested]^ Since 
such a magnetic transition in this temperature region 
is not refiected in the results of either neutron scatter- 
ing or susceptibility experiments ^i^ii it may be assumed 



that the inferred magnetic order does not involve long- 
range order and static internal fields. Considering the 
electronic configuration of the Pr'^"*" ions, the above re- 
sult is really unexpected and calls for additional checks. 
Since the ^SR results may be due to dynamically-induced 
non vanishing internal fields, it is important to use a 
microscopic technique for investigating the phenomenon 
further. In this paper we report a systematic study of 
PrCu2 by means of Cu-based nuclear magnetic resonance 
(NMR) and nuclear quadrupole resonance (NQR) exper- 
iments between 4.2 K and room temperature. 



II. EXPERIMENTAL 

Nuclear resonance measurements were made us- 
ing standard spin-echo techniques employing a phase- 
coherent pulsed spectrometer. ^^^^^Cu NMR and NQR 
spectra were collected by measuring the integrated spin- 
echo signal as a function of the exciting radio-frequency 
(RF) in fixed external magnetic fields. The spin-lattice 
relaxation rate (SLRR) T]~^ was measured by destroying 
the z-component of the equilibrium nuclear magnetiza- 
tion by means of a short RF pulse and monitoring its 
recovery as a function of a variable delay between this 
preparation pulse and a spin-echo sequence. The spin- 
spin relaxation rate (SSRR) T2~^ was obtained from the 
T dependence of the signal after a 7r/2 — r — tt sequence. 

Single crystalline PrCu2 samples were grown at Os- 
aka University by means of the procedure described in 
Ref. [nl. Since the system is metallic, the NMR signal 
from the single crystal is very small because the penetra- 
tion depth of the RF pulses is of the order of 1 ^m only. 
We observed an unexpected splitting and broadening of 
the NMR fines, which we traced back to a non-perfect 
orientation of the applied magnetic field with respect to 
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Figure 1: ^^'''^Cu-NMR powder spectra of PrCu2 at selected 
temperatures and in a field of 7.049 T. The vertical dashed 
lines mark the reference NMR frequencies for the two Cu- 
isotopes.— 



one of the crystalline axes of the sample. From a de- 
tailed analysis we concluded that meaningful NMR mea- 
surements on a single crystal of PrCu2 require a preci- 
sion in sample orientation that is much better than ±1°, 
which could not be achieved with our experimental set- 
up. For this reason we powdered the single-crystalline 
sample which not only served to avoid the orientation 
problem but also helped to enhance the NMR signal. 



III. NMR MEASUREMENTS 

Figure [T] shows the Cu-NMR spectra of PrCu2 pow- 
der at selected temperatures. The 300 K spectrum is 
characterized by two main features centered at 80 and 
86 MHz, which are ascribed to the NMR central-line of 
the ^"^Cu and ^^Cu isotopes, respectively. Both isotopes 
possess a 3/2 nuclear spin but the gyromagnetic ratio of 
^^Cu is roughly 7% larger than that of ''^Cu. We ob- 
serve no sign of any contribution from the ^^^Pr nuclei 
(reference frequency ^ 92 MHz at 7.049 T), most prob- 
ably because their extremely large hyperfine coupling^ 
results in very short relaxation times. Both Cu NMR 
lines are split by ~ 1 MHz. We ascribe this splitting 



to a second-order quadrupolar perturbation of the Zee- 
man Hamiltonian. The wings observed between 75 and 
90 MHz correspond to the first-order quadrupolar per- 
turbation. These wings are more prominent at lower 
temperatures and it may be noted that, rather than a 
sharp edge, they exhibit a smooth edge suggesting, as 
explained below, a non-zero anisotropy ry of the electric 
field gradient (EFG). With decreasing temperatures, the 
complexity of the spectrum increases. As may be seen in 
fig. [H the asymmetry and width of the lines grows and 
they partially merge together. The quadrupolar wings 
are strongly affected as well and a large overall shift of the 
lines towards higher frequencies is observed, indicating a 
progressively increasing Knight-shift (K). The strongly 
asymmetric line-shape suggests a large K-anisotropy. 

The very complex structure of the low temperature 
NMR spectra, caused by the simultaneous EFG and K 
anisotropics, makes a quantitative analysis of the spec- 
tra very difficult. Therefore, we keep the interpretation 
mostly at a qualitative level. To this end we performed 
several simulations of the ^'^'^^Cu NMR powder-spectrum 
for different values of the microscopic parameters char- 
acterizing K and EFG. The simulations involve a full di- 
agonalization of the nuclear Hamiltonian comprising the 
sum of the quadrupolar and the Zeeman term corrected 
by K. The spectrum is then obtained from the calcu- 
lated transition energies and intensities, averaged over 
random crystal orientations in order to simulate a pow- 
der spectrum. The microscopic parameters characteriz- 
ing the quadrupolar Hamiltonian are the largest com- 
ponent Vzz of the EFG and the anisotropy 77 in the xy 
plane. The parameters describing K are the isotropic and 
anisotropic components i^iso and i^ani, respectively, as 
well as an axial-asymmetry term e. The notation is stan- 
dard and exemplified, for instance, in Ref. [13 (pages 62- 
64). Figure [2lja) shows the simulated spectra that we ob- 
tained by considering only an isotropic i^igo ~ 1.5% and 
no quadrupolar contribution. Naturally, this spectrum 
consists of only one line for each isotope, slightly shifted 
with respect to the reference frequency. Figure [2][b) was 
obtained by adding an EFG with Vzz = 0.32 a.u. and 
?7 = 0. It is clear that this term induces the appear- 
ance of the wing structures and a splitting of the central 
lines. The effect of adding a non-zero anisotropy 77 is 
demonstrated in figs.[2ljc) and (d) where the calculations 
with 77 = 0.10 and 77 = 0.20, respectively, are displayed. 
Upon enhancing 77, the sharp edge of the wings is pro- 
gressively smeared out, leading to a simulated spectrum 
with a shape resembling the experimental data at 300 K. 
The overall comparison of the experimental data with the 
simulations of fig. [2ljb-d) suggests that the actual compo- 
nents of the EFG acting at the copper site in PrCu2 are of 
the order of Vzz = 0.32 a.u. and 77 = 0.10, at least down 
to 100 K. An ah initio calculation in which the internal 
coordinates of the 12 ions in the unit cell of PrCu2 were 
optimized assuming the experimental lattice parameters 
(a = 4.4087 A, 6 = 7.0551 A, and c = 7.4441 kf- using 
the DMol^ code^ yielded Vzz = 0.48 a.u. and -q = 0.11. 
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Figure 2: Simulated ' Cu-NMR powder spectra in a field 
of 7.049 T with Ki^o = 1-5% and for different values of the 
EFG parameters (see text). 



Figure 3: Simulated powder ' Cu-NMR spectra at a field 
of 7.049 T and for different values of the K parameters (see 
text). 



Although the calculated absolute value of the EFG is 50% 
larger than that obtained by comparing experiment and 
simulations as described above, the direction of the dom- 
inant component of the EFG tensor is clearly established 
to lie in the be plane. This is important for the discussion 
in section ITVB I 

Figure [3ljb) shows a simulated spectrum where the 
isotropic shift of 1.5% is augmented by an anisotropic 
component i^ani = 1.6% but no axial-asymmetry. The 
non-zero i^ani induces an asymmetric broadening of the 
central lines keeping rather sharp edges. The e value is re- 
sponsible for the position of the maximum with respect 
to the edges, as shown in figs. [HJc) and (d), where the 
simulated spectra with e = 0.6 and e = 2.0 are shown, 
respectively. 

Although these simulations allow for a qualitative de- 
scription of the main features of the NMR spectra, the 
situation turns out to be much more complicated when 
we try to simultaneously take into account the influence 
of both EFG and K. These two contributions are indeed 
both non-negligible and, in principle, temperature de- 
pendent. Moreover the corresponding tensors are usu- 
ally non diagonal in the crystalline-axes frame and their 
main frames, i.e. the frames in which the two tensors 
are diagonal, are not necessarily oriented along the same 
directions. The local symmetry of the Cu-sites provides 
one constraint for the EFG and K main frames. They 
should have one principal axis parallel to a. Concern- 
ing the orientation of the other two axes in the be plane, 
nothing can be concluded from considering the crystal- 
structure. Therefore the angle between the K and EFG 
main frames in the plane is unknown. All the above con- 
siderations lead to the conclusion that the influences of 
both K and EFG on the spectral features cannot be re- 
garded independently. 



In spite of the difficulties described above, it is still 
possible to extract some quantitative information from 
the NMR spectra shown in fig. [TJ First of all, we can 
estimate the quadrupolar frequency of ^^Cu from the 
widths of the wings. The data imply =8 — 9 MHz 
with no signiflcant change upon decreasing temperature, 
at least down to 50 K. Inserting Vzz = 0.32 a.u. and 
?/ = 0.10 into the relation 



ieVzzQ 
2hl{21 - 1 



(1) 



where Q is the nuclear quadrupole moment and / is the 
nuclear spin, we obtain i^g = 8.3 ±0.8 MHz. The second 
quantity that we can estimate is the largest principal 
component /•sTmax of the K tensor. In the presence of 
a K anisotropy, the spectrum adopts the typical shape 
displayed in flgs.[3Kb)-(d).— The edges of the asymmetric 
broadening occur at frequencies (1 -I- Kyain)'^o and (1 -I- 
-?^max)t'0; where i^min is the smallest component of the K 
tensor and vq is the reference frequency of the considered 
nucleus.— The position of the resonance maximum is at 
(1-1- i^int)'^0; where Kj^t is the third (intermediate) K- 
tensor main-component.— The temperature dependence 
of the high-frequency edge can be followed down to ~ 
20 K, thus allowing us to estimate -ftTmaxlT")- 

We assume that the hyperflne fleld Bq^ acting at the 
Cu-site is mainly induced by the polarization of the Pr 
magnetic-moments Jlp^ via an interaction mediated by 
the conduction electrons, such that 



p>hf 

^Cu 



(2) 



where A is the tensor associated with the transferred hy- 
perflne fleld. The orientation of one principal axis of the 
K tensor along the a axis and the fact that it represents 
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the easy magnetization direction^ indicate that the max- 
imum K component most likely lies along a, as was also 
suggested by our preliminary measurements on the single 
crystal (data not shown here) . The same is implied by our 
estimate of K^^axiT), which is found to be proportional 
to the susceptibility Xa{T) along the a axis. For a given 
external magnetic field Bq, and /Jpr can directly be 
obtained from K and the susceptibility, respectively. Us- 
ing eq. [2] we can calculate the diagonal term of A along 
a, namely 

Aaa^KaiT)/Xa{T), (3) 

which is, of course, temperature independent and has a 
value of 0.54 ± 0.02 T//Ub. This value is very useful in 
the discussion of the NQR data. 

A temperature-independent Aaa and the fact that we 
do not observe significant changes of the NMR spectrum 
around 65 K are not compatible with an ordering of the 
Pr-moments at that temperature. Because of the inter- 
action between the Pr-moments and the conduction elec- 
trons, any alignment of the Pr moments is expected to 
have a strong influence on the Cu NMR response. 



IV. NQR MEASUREMENTS 
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A. Spectra 

The rather large value of vq indicated by the NMR 
data suggests that an NQR signal can be observed in 
an accessible frequency range. Figure [4] shows Cu NQR 
spectra at selected temperatures. The recorded data re- 
flect the expected simple NQR spectrum with one single 
line per isotope, reflecting the respective ±1/2 — > ±3/2 
transitions. The relative peak positions of the ^^Cu and 
^^Cu signals are consistent with the nuclear quadrupole 
moment, which is roughly 8% larger for ^^Cu than for 
^^Cu. From the raw data it is clear that the '^^Cu NQR 
frequency is within 10% of the estimate of vq obtained 
from the NMR spectra. The NQR frequency increases 
with decreasing temperature, whereas no major changes 
are observed in the T-dependence of the line-widths. In 
order to obtain more quantitative information, we ap- 
plied a best-flt procedure that is based on the frequency 
dependence of the intensity S{i') given by 



= 



Aa 



65 



Aaes 



(r/2)2 



(z.-i.Q)2 + (r/2)2 



(4) 



where A is the overall intensity, ags^es are the natural 
abundances of ^^■^'^Cu, uq is the ®^Cu NQR frequency, 
r is the line-width (assumed to be equal for both iso- 
topes), and (3 = 1.078 is the ratio between the quadrupo- 
lar moments of ^"^Cu and ^^Cu. The experimental inten- 



Figure 4: 



'Cu-NQR spectra of PrCu2 at selected temper- 



atures (symbols). Solid lines are best-fit curves to the data 
(see text). 



sity I^v, T) is 

I{v,T) = S{v)u 



1 — exp 



hv 



(5) 



where the multiplicative factor v is due to the fact that 
our measurement is inductive; the term in square brack- 
ets represents the Boltzmann population factor. With 
this procedure, we fltted all the spectra using only A, 
Vq, and r as free parameters. Physically relevant are the 
latter two and we plot them as a function of temperature 
in flg. [5l 

The frequency vq increases linearly with decreasing 
temperature down to Tjt , where a sudden drop in vq (T) 
is observed. This reflects the JT structural transition 
which induces a sudden change in the EFG and thus in 
i^Q. The overall change of vq between 300 K and 20 K is 
about 5%. Since vq is related to both Vzz and rj by means 
of eq.[Tl it is not possible to sort out whether the observed 
change in vq is to be ascribed to Vzz, or 77, or both. Nev- 
ertheless, if we assume a constant Vzz, the vq{T) values 
imply that the anisotropy should change from 77 ^ 0.10 
at 300 K to approximately 0.60 at 20 K. Such a dra- 
matic change in the asymmetry is very unlikely because 
it would require a strong modiflcation of the electronic 
density, which does not seem to be compatible with the 
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Figure 5: Symbols: (a)'^^Cu-NQR frequency and (b) line- 
width of PrCu2 as a function of temperature, as obtained from 
the best-fits to the NQR spectra. The dashed lines represent 
(a) a linear fit and (b) the average of the data above 7 K. 



small anisotropy of the thermal expansion coefficients. — 
Moreover no evidence of such a large variation can be in- 
ferred from the NMR data shown in fig-HJ because in the 
presence of a large rj the edges of the quadrupolar wings 
are dramatically smeared out. It is therefore reasonable 
to ascribe most of the change in i^q to a corresponding 
increase of Vzz , probably induced by the thermal contrac- 
tion of the lattice.— The lack of any anomaly in vq{T) 
around 65 K does not add support to the magnetic-order 
scenario inferred from the ^SR data. 

Between room-temperature and Tjt the line-width T 
is almost constant. In particular, also r(T) exhibits no 
anomaly around 65 K. Likewise the absolute values of T 
are in contrast with a magnetic order of the Pr moments. 
Because of the interaction of the Pr-moments with the 
conduction electrons, magnetic order among them would 
imply a local field acting on the Cu-sites. As dis- 
cussed above (sect. IIIip . the contribution to this field 
of any Pr moment oriented along a can be calculated 
by means of the corresponding hyperfine constant, i.e. 
-^Cu ~ ^aaMPr- We Calculated the local magnetic fields 
acting on Cu sites by invoking the six nearest-neighbors 
Pr moments, adopting the claimed AF configuration with 
an incommensurate modulation of the moments. The 
corresponding vector ~q = (0.24,0,0.68) and the modu- 



lation amplitude is 0.29 /xbi^ We obtain a distribution 
of B^^^ ranging from to 0.06 T. Such a distribution of 
fields should be refiected in a broadening Av = 1.4 MHz 
of the NQR line. This value is almost one order of mag- 
nitude larger than T of our signals, thus indicating that 
a magnetic order of the claimed type^^ is incompatible 
with our data. 

Since the features extracted from the NQR spectra do 
not confirm, but rather rule out AF ordering of the Pr 
moments at 65 K, it seems in order to discuss possible 
causes for this discrepancy. At present we consider two 
scenarios which we discuss in more detail in section [Vl 
The first is based on the assumption that the apparent 
order observed with muons is mainly due to a muon- 
induced local enhancement of slow spin fiuctuations. As 
a second possibility we consider that the observed order 
is not static, but fiuctuates on a time-scale of the order 
of 1 /IS. With /iSR and the related precession frequency 
of the order of 10 MHz,— a non-zero field would be in- 
dicated. The characteristic time period in employing the 
NMR/NQR spin-echo technique are in the range of tens 
of microseconds and therefore the field fiuctuations in the 
MHz range would be averaged to zero and thus would not 
affect the resonance signal. Under these circumstances, 
/iSR would indicate a "static" magnetic correlation with 
a distribution of internal fields, whereas in NMR/NQR 
the fields would be averaged to zero. Such a slow dynam- 
ics has no simple explanation. Of course, the presence of 
a temporally fiuctuating field should be refiected in the 
dynamical magnetic properties on a microscopic scale, 
namely the SLRR Tf ^ and the SSRR Tj"^- At any rate, 
the onset of magnetic-order is, usually in the vicinity of 
the transition, accompanied by strong magnetic fiuctua- 
tions, which should be refiected in the relaxation rates. 



B. Relaxation 

With the procedures described in section [Ul we mea- 
sured the time recovery of the longitudinal mL{t) and 
transversal mT(t) magnetization. Representative exam- 
ples of raL{t) and ■mT{t) are shown in the insets of 
fig-EJa) and (b), respectively. All the mL[t) curves can 
be fitted by means of the standard exponential function 
that applies for 1 = ^/2^ namely 



mL{t) = mL(oo) 



1 -exp ( - — 



(6) 



where 771^(00) is the equilibrium magnetization, Ti is the 
spin-lattice relaxation time, and p is an 77-dependent pa- 
rameter, which can be approximated with p ~ 3 for < 
0.2.— Because the standard dipole-dipole interaction 



V 



< 



is not the main spin-spin relaxation mechanism, mT(t) 
cannot be fitted by means of a simple exponential func- 
tion. The stronger indirect interaction mediated by the 
conduction electrons via the Ruderman-Kittel-Kasuya- 
Yosida mechanism leads to a gaussian correctioni^ii^ 
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Temperature (K) 

Figure 6: (a) Spin-lattice and (b) spin-spin relaxation rates 
as a function of temperature. Dashed lines in (a) are power- 
law fits to the data in the temperature range above and below 
100 K. The uncertainty of the data points in (a) is of the order 
of the diameter of the open circles. The horizontal dashed 
line in (b) is a guide to the eyes. Insets: examples of the time 
dependent longitudinal (a) and transversal (b) magnetization. 
Solid lines in the insets of (a) and (b) are best fits to data 
according to equations [6] and respectively. 



such that 



mT{t) ~ tot(O) exp 



t 



2 \T2gJ 



(7) 



where tot(O) is the initial transverse magnetization, T2G 
is the gaussian spin-spin relaxation time, and Tr = 
3Ti/(2 -I- r) is a time-scale depending on the SLRR and 
on its anisotropy r. Since this anisotropy is expected 
to be not as large as, e.g., in the layered copper oxides 
(r ^ 4) , we assume r ~ 1 and hence T/j ^ Ti . As long as 
Tr is significantly larger than T2G, their values are quite 
independent. Our assumption = Ti is thus justified 
a posteriori, both by the good quality of the fits and by 
the fact that we obtain systematically T2G < Ti. 

Prom the fits based on equations [6] and [71 we obtain 
the temperature dependencies of the SLRR and SSRR, as 
shown in figs.[6lja) and (b), respectively. Neither T^^{T) 
nor T2g{T) exhibit the expected features that would re- 
fiect an onset of magnetic-order around 65 K. The fiuc- 
tuations that accompany such a transition should cause 



a distinct enhancement of the relaxation rates, result- 
ing in a peak in T^^{T) around the ordering tempera- 
ture. From fig. m^b) we take it that T^QiT) is constant 
down to around 40 K where it starts to increase with 
decreasing temperature, displaying a peak at Tjt- The 
decrease of T^q (T) below Tjt most Hkely refiects the re- 
duction of the quadrupolar fiuctuations. Also Tj~^(T) 
peaks at Tjt, but the behavior at higher temperature is 
more difficult to interpret. The temperature dependence 
of the SLRR is neither compatible with a dominating re- 
laxation via conduction electrons {T^^{T) T), nor a 
relaxation via local magnetic moments {T^^ = const.) 
or by spin- wave type excitations {T^^ ~ T^). Instead 
Tf ^ (T) exhibits a power-law type behavior characterized 
by a negative exponent (-1.4) between 8 and 100 K and 
a positive one (0.4) above 100 K. In order to interpret 
this puzzling behavior, we tried to determine the leading 
spin-lattice relaxation mechanism by measuring Ti also 
for 65 Cu. With ^^^^Ti, 63.65^, and ^^.esg denote the 
spin-lattice relaxation time, the gyromagnetic ratio, and 
the quadrupole moment of 63'65(2;y^ respectively. If the 
dominant mechanism is magnetic, 



63r 



Ti/i^^Ti = r-fT-'iy = 1.167, 



(8) 



whereas for a relaxation driven by the quadrupolar inter- 
action we expect 



63r 



Ti/^^Ti = C^Qr^QY = 0.872. 



(9) 



We estabhshed that below 100 K esy^/esj^^ ^ n 
indicates that the relaxation is predominantly of mag- 
netic origin. At higher temperature we find ^^Ti/^^Ti ~ 
1, suggesting that the magnetic and quadrupolar inter- 
actions are of comparable importance. The fact that the 
magnetic fiuctuations dominate below 100 K again sug- 
gests that the internal field, if any is present, is not static. 
It is also worth mentioning that the puzzHng T^^{T) 
cannot be ascribed to the temperature dependence of 
the population of the Pr levels in simple terms. Such 
a case would require that the lifetimes of these levels are 
different for each of them. Bearing in mind the single- 
exponential behavior exhibited by all our mL{t) curves, 
this possibility appears as rather unlikely. 



V. COMPARISON BETWEEN NMR/NQR AND 
/iSR MEASUREMENTS 

Our ab initio calculations indicate that the dominant 
axis of the EFG and thus the quantization direction in 
NQR, is orthogonal to the a axis (see section HlH) . If we 
assume that the internal field along a observed by ^SR 
is not static but fiuctuates, this would result in a large 
SLRR whereas it should have only minor effects on the 
SSRR. This is consistent with our observation that on de- 
creasing the temperature, T^^{T) increases much more 
rapidly than T2q{T). As already mentioned (sect. lIV A|) . 
this scenario could reconcile the NMR/NQR data with 
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Figure 7: Sketch of the positions of the Pr ions surrounding 
the muon, as obtained from the ab initio calculations (see 
text). The arrows indicate the displacement of the Pr ions 
(multiplied by 50) induced by the presence of the muon. 



site which is a factor 6.6 larger than the average one. Re- 
ducing the imposed spin polarization by a factor 4 yields 
a spin density at the muon 12 times larger than the aver- 
age value. One might argue that the conclusions drawn 
from a calculation in which a hydrogen atom is placed in 
every second cell cannot be appHed to the analysis of /xSR 
experiments, which are performed with isolated muons. 
Our experience shows that if we double the number of 
hydrogens to one per cell, the enhancement of the spin 
density at the muon-site is reduced by 15%. In other 
words, a larger dilution is expected to enhance the influ- 
ence of the muon. So there is no doubt that the presence 
of the muon affects the response of the system, namely a 
small spin polarization in the system will be amplified at 
the muon site. A quantitative assessment of the temper- 
ature dependence of this extrinsic effect is out of reach 
of the presently used calculation method. 



the /^SR results, provided that the time-scale of the spin 
polarization is long 1 ^s). This would also be consis- 
tent with the observation of a gaussian transverse relax- 
ation, since the presence of such a term requires an en- 
hancement of the low-frequency components of the local 
susceptibility,— i.e., a slowing down of the spin dynamics. 

In order to reconcile the results of the NMR/NQR and 
/LtSR experiments, we also consider the possibility that 
the magnetic order observed by ^SR is mimicked by a 
quasi-static local enhancement of moment polarization 
by the probe itself. The muon, considered as a posi- 
tive point charge, will attract electrons and thus lead to 
a rearrangement of charge and possibly of the ion po- 
sitions. If, because of a slow enough spin fluctuation, 
the intrinsic electronic charge distribution is spin polar- 
ized during the duration of the muon lifetime, the muon 
will experience a spin density exceeding the one induced 
by the fluctuation. For our calculationsi^ we introduce 
a hydrogen atom^fi at the position of the muon derived 
by Schenck et al. (1/2,1/4,0.6406)^ in every second 
cell along the a direction, which amounts to one muon 
per 8 formula units, and let all atoms relax at flxed vol- 
ume. The resulting coordinates for the position of the 
muon are (1/2, 1/4, 0.6465), independent of the imposed 
spin polarization between and 1 per unit cell. The 
Pr-ions in the nearest neighbor triangle are displaced as 
shown in flg.[7l The Pr-Pr distance is reduced by 0.3% 
which corresponds to a compression by approximately 
5 kbari^ At this point, it is worth mentioning that at 
external pressures exceeding 12 kbar, AF ordering of the 
Pr moments has been observed for T < 9 K,i2L2^ thus 
indicating that a lattice compression enhances the mag- 
netic correlations. Imposing a spin polarization of 1 hb 
per unit cell, which corresponds to an average spin den- 
sity of 6.4 • 10"** a.u., leads to a spin density at the muon 



VI. CONCLUSIONS 

The results of our NMR/NQR study of both the static 
(spectra) and dynamical (relaxation rates) microscopic 
magnetic properties of PrCu2 provide no evidence for the 
onset of any magnetic order around 65 K that was pre- 
viously claimed from ^SR dataji^ As argued above, the 
/iSR results may not reflect magnetic order but are rather 
due to a muon-induced local enhancement of slow polar- 
ization fluctuations. If the fluctuation has a time scale of 
the order of 1 /xs, it could, as argued in section HV A^ ap- 
pear as static in ^SR and vanishing in NMR/NQR. The 
origin of the anomalous power-law behavior of the SLRR 
is not yet understood, but it clearly indicates that the 
magnetic behavior of PrCu2 is far from that of a simple 
Pauli, Curie- Weiss type, or spin-wave system. 

The present results call for further studies in order to 
clarify the complex magnetic properties of PrCu2. In 
particular it will be important to carry out NMR/NQR 
measurements in an AF phase. This could be obtained 
either by means of measurements below 54 mK or at 
pressures exceeding 12 kbar, where AF ordering among 
localized Pr moments has been observed for T <9 
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